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Abstract 

Device-to-device (D2D) communication enables direct communication between proximate devices 
thereby improving the overall spectrum utilization and offloading traffic from cellular networks. This 
paper develops a new spatial model for D2D networks in which the device locations are modeled as 
a Poisson cluster process. Using this model, we study the performance of a typical D2D receiver in 
terms of coverage probability under two realistic content availability setups: (i) content of interest for a 
typical device is available at a device chosen uniformly at random from the same cluster, which we term 
uniform content availability, and (ii) content of interest is available at the closest device from the 
typical device inside the same cluster, which we term k-closest content availability. Using these coverage 
probability results, we also characterize the area spectral efficiency (ASE) of the whole network for the 
two setups. A key intermediate step in this analysis is the derivation of the distributions of distances 
from a typical device to both the intra- and inter-cluster devices. Our analysis reveals that an optimum 
number of D2D transmitters must be simultaneously activated per cluster in order to maximize ASE. 

This can be interpreted as the classical tradeoff between more aggressive frequency reuse and higher 
interference power. The optimum number of simultaneously transmitting devices and the resulting ASE 
increase as the content is made available closer to the receivers. Our analysis also quantihes the best 
and worst case performance of clustered D2D networks both in terms of coverage and ASE. 
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I. Introduction 

Enabling direct communication between devices located in close proximity, termed device-to- 
device (D2D) communications, has several benefits compared to the conventional approach of 
communicating through a base station in a cellular network Q-Q. First, the spectral efficiency 
of the direct link is typically much higher due to a smaller link distance. Second, this circumvents 
the need to establish an end-to-end link through a base station, thereby offloading traffic from 
cellular networks. Third, while the D2D network can be visualized as an ad hoc network, it 
incurs a much lower protocol overhead due to the assistance it gets from the existing cellular 
network. All these benefits make it an attractive component for both the current 4G and the 
future 5G networks 0-0. Clearly, content centric nature of D2D communication opens up 
several exciting possibilities that were not quite possible with traditional cellular architecture. 
This is primarily driven by the spatiotemporal correlation in the content demand @-[10|. In 
particular, when a device downloads a popular file, it can deliver it locally to its proximate 
devices whenever they need them pT|-pA|. We term each such set of proximate devices as a 
cluster. The performance of a typical D2D link within each cluster will mainly depend upon 
where the content of interest is available with respect to the typical receiver and the number of 
other D2D links active in the network. Comprehensive modeling and performance analysis of 
this clustered D2D network using tools from stochastic geometry is the main goal of this paper. 


A. Motivation and Related Work 

Modeling and analysis of D2D communication has taken two main directions in the literature. 
The first one focuses on characterizing the scaling of per-device throughput as a function of 


the network size; see [12|-[14| for a small subset. To maintain analytical tractability, protocol 
model is typically assumed under which the transmission between two devices is successful 
only if (i) the distance between them is smaller than a certain predefined value, and (ii) there 
are no other active transmitters in the immediate neighborhood of the receiver. While these 
assumptions are somewhat restrictive, more general results are typically obtained by means of 
extensive simulations. The second direction, which is also more relevant to our work, focuses on 
characterizing metrics, such as the distribution of per-device throughput and coverage probability, 
using tools from stochastic geometry under more general physical layer models in which the 
metrics are defined in terms of the actual received powers from the desired and interfering 
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devices, as opposed to Euclidean distances that appear in the protocol model discussed above. 
While stochastic geometry has successfully been applied to study various aspects of ad hoc and 
cellular networks over the past decade; see p3|-p0| for a small subset, it has also been used 


more recently to study D2D networks pT|-[29|. We discuss these works in more detail next. 
Mode selection in D2D-enabled uplink cellular networks where traffic can be offloaded from 


the cellular network to the D2D network is studied in [21|-[23|. Interference management 


between the cellular and D2D networks is investigated in [24|-[28|. The performance of the 
multicast transmission wherein each D2D transmitter (D2D-Tx) has a common massage for 


all the intended D2D receivers (D2D-Rxs) inside the cluster is analyzed in [29|. The common 
approach in all these works is to model the locations of the the D2D-Txs as a Poisson Point 
Process (PPP) while two approaches are considered for modeling the locations of the D2D- 
Rxs. In the first approach, to lend analytical tractability, the network is modeled using a Poisson 


Dipole Process (PDP) where the D2D-Rxs are located at a fixed distance from the D2D-Txs pi |, 
p4|-p^. Although this is a good first order model, the assumption of fixed link distance is 
quite restrictive. This assumption is relaxed by assuming that the intended D2D-Rx is uniformly 


distributed within a circle around its serving D2D-Tx [22|, [23|, [29|. However, neither of 
these stochastic geometry-based approaches captures the possibility of having multiple proximate 
devices any of which can act as a serving device for a given device, which is quite fundamental to 
D2D networks pT|-[ 131. In this paper, we address these shortcomings by developing a new and 
more realistic spatial model for D2D networks in which the devices form clusters. We consider 
out-band D2D in which D2D and cellular transmissions do not interfere with each other. More 
details of the model along with the other main contributions of this paper are provided next. 


B. Contributions and Outcomes 

Realistic tractable model for D2D networks. We develop a new and more realistic way of 
modeling D2D networks in which the device locations are modeled as a Poisson cluster process. 


in particular a variant of a Thomas cluster process |15|. This is unlike the popular approaches 
where the device locations are assumed to be uniform over the plane, such as in the PPP and 
PDP models discussed above. The proposed model captures the fact that a given device typically 
has multiple proximate devices any of which can potentially act as a serving device. Using tools 
from stochastic geometry, we characterize the performance of this D2D network for two content 
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Fig. 1: Proposed D2D cluster model where devices are normally distributed around each cluster center. 


availability scenarios: (i) uniform content availability, where the eontent of interest for a typieal 
deviee is available at a deviee ehosen uniformly at random from the same eluster, and (ii) k- 
closest content availability, where the eontent of interest is available at the elosest deviee to 
a typieal deviee inside the same eluster. It should be noted that while Poisson eluster proeess has 


been used in the literature to model wireless networks, in partieular see [30|, the performanee 
analysis is usually performed at a point that may not necessarily be a part of the point process. 
Unlike this seenario, we perform analysis at a typical device, whieh by definition is a part of the 
cluster proeess. This setup brings forth new teehnical ehallenges, e.g., the need to eharaeterize 
the distribution of distances from the typieal deviee to other deviees, as discussed below. 

New distance distributions enabling ASE and coverage probability analysis. We derive exaet 
expressions for coverage probability of a typical device and ASE of the whole network for 
the two scenarios described above. Several elosed-form bounds and approximations are also 
derived. As intermediate results, we eharaeterize the distributions of the distanees of the typieal 
device to its serving device and intra- and inter-cluster interfering deviees in the two cases. 
The key enabler of our analysis is the observation that conditioned on the appropriate random 
variables, the distanees from a typieal deviee to any other randomly ehosen deviee in a Thomas 
cluster proeess are independent and Rician distributed. The fc-closest content availability case 
additionally requires the use of order statistics in order to characterize the serving distanee 
distribution and a more eareful treatment of the intra-eluster interferenee field. Using these 
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distributions, we derive Laplaee transforms of both intra- and inter-eluster interferenee powers, 
using whieh the eoverage probability and ASE are eharaeterized. 

System design insights. Our analysis leads to several system design guidelines. First, it reveals 
the existenee of the optimal number of links that must be aetivated per eluster in order to max¬ 
imize ASE. This ean be interpreted as the elassieal tradeoff between more aggressive frequeney 
reuse and higher interferenee power. The optimum number of simultaneously aetive links per 
eluster and the resulting ASE inerease as the eontent is made available eloser to the reeeivers. 
For typieal operational regimes of interest for D2D networks, our results reveal that signifieant 
gains ean be aehieved by aetivating optimum number of links eompared to strietly orthogonal 
strategy in whieh only one link per eluster is aetive. The elosest eontent availability strategy 
also allows to eharaeterize the best and worst ease performanees of a elustered D2D network in 
terms of eoverage probability and ASE by tuning the value of k, where the best and worst eases 
eorrespond to the when the eontent is available at the elosest and farthest deviee, respeetively. 

II. System Model 

We assume that eaeh deviee has a eertain eontent that ean be requested by the other deviees 
in the same eluster and the deviees aeross elusters need not oommunieate. This ean be justified 
in a praetieal setup from multiple perspeetives, two of whieh are: (i) the inter-eluster distanees 
will be typieally mueh larger than the intra-eluster distanees between deviees, thus making it 
easier to oommunieate within a given eluster, and (ii) the deviees in one eluster may not have 
information that is of interest to the deviees in the other elusters. For example, deviees forming 
a eluster in a sports bar are more likely to be interested in sports-related eontent as opposed to 
deviees in an aeademio setting, suoh as a library, where the nature of popular eontent may be 
entirely different. While the tools developed in this paper ean be extended to handle the ease 
where oommunioation aeross elusters is allowed, it is not in the seope of this paper and is left 
as a promising future work. We now diseuss the key modeling details for this setup. 

A. Spatial Setup and Key Assumptions 

The loeations of the deviees are modeled by a Poisson cluster process, where the parent point 
proeess is modeled by a PPP <I)c with density Ac, and the offspring point proeesses (one per 


parent) are eonditionally independent [311. The union of all the offspring points eonstitutes a 
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Poisson cluster process. The parent and offspring points will be henceforth referred to as the 
cluster centers and the cluster members (or simply devices), respectively. The cluster members 
are assumed to be independent and identically distributed (i.i.d.) according to a symmetric normal 
distribution with variance cr^ around each cluster center x G $c- Therefore, the density function 
of the device location ?/ G relative to a cluster center is 


fvidj) 


1 


27r cr^ 


exp 



( 1 ) 


If the number of devices in each offspring process, i.e., cluster, were Poisson distributed, this 


process is simply a Thomas cluster process p3| . However, to simplify certain order statistics 
arguments in the sequel, we assume that the total number of devices per cluster is fixed and 
equal to N. As will be evident from the discussion below, the number of simultaneously active 
devices will still be different across clusters, thereby providing sufficient generality to the model. 

The proposed model is illustrated in Fig.[^ The set of all devices in a cluster a; G <I>c, denoted 
by is partitioned randomly into two subsets: (i) set of possible transmitting devices denoted 
by and (ii) set of possible receiving devices denoted by The set of simultaneously 
transmitting devices in this cluster is denoted by C A/f, where \B^\ is assumed to be Poisson 
distributed with mean fh conditioned on \B^\ < Note that in the limiting case, about half of 
the devices in each cluster will transmit to the other half. Therefore, for notational convenience, 
we assume that the total number of transmitting devices per cluster is limited to M = N/2. 
This can be easily relaxed in case one wants to allow all N devices to be transmitting in a 
certain application. Without loss of generality, we perform analysis for a typical device, which 
is a randomly chosen device in a randomly chosen cluster, termed representative cluster, inside 
the network. Assuming the cluster center of the representative cluster to be located at xq G <I>c> 
the typical device by definition is in . Due to the stationarity of this process, we assume that 
the typical device is located at the origin. Since the performance of the D2D link to this receiver 
depends upon where in the cluster is the data available, we consider following two setups to 
model content availability in the representative cluster: 


1) Uniform content availability. The content of interest for the typical device is available at 
a device that is chosen uniformly at random from the set in the same cluster. 

2) k-closest content availability. The content of interest to the typical device is located at 
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its closest device from the set in the same cluster. By tuning the value of k, the 
content can be biased to lie closer (small k) or farther (large k) from the typical device. 
After fixing the location of the serving device as per one of the strategies above, the intra-cluster 
interfering devices are sampled uniformly at random from the remaining M — 1 devices in 
in the representative cluster. Since a representative cluster has a serving device by definition, 
for concreteness we assume that the number of interfering devices is Poisson distributed with 
mean m — 1, which means that a total of fh devices are active on average in this cluster. 
Similarly, the inter-cluster interfering devices are sampled uniformly at random from the set of 
transmitting devices of each cluster, such that the number of active devices in each cluster is 
Poisson distributed with mean fh conditioned on the total being less than M. 

B. Channel Model 

Recall that the cluster center of the representative cluster is assumed to be located at xq G <fic- 
Assume that D2D-Tx of interest (serving transmitter) is located at yo inside that cluster. The 
distance of this D2D-Tx from the typical device (D2D-Rx of interest) is denoted by r = |||/o+a^o||> 
where r is a realization of random variable R whose distribution depends upon the content 
availability strategy as discussed in detail in the sequel. Assuming transmit power of each device 
to be Pd, the received power at D2D-Rx of interest is 

P = Pd/io|ko + l/oir“, (2) 

where ho ~ exp(l) is i.i.d. exponential random variable which models Rayleigh fading and a is 
path loss exponent. Incorporating shadowing is left as a promising direction of future work. To 
define interference field, it is useful to define the set of all simultaneously active D2D-Txs as: 

'km = (3) 

where recall that is the set of simultaneously active D2D-Txs inside a cluster x E <kc. In this 
network, the total interference caused at the D2D-Rx of interest can be written as the sum of two 
independent terms: (i) intra-cluster interference caused by the interfering D2D-Txs inside the 
representative cluster, and (ii) inter-cluster interference caused by simultaneously active D2D- 
Txs outside the representative cluster. Recalling that the D2D-Tx of interest is located at yo with 



respect to the cluster center xq, the intra-cluster interference power can be expressed as: 


-^Tx—cluster ^ ^ “f VW • (4) 

yeB^Q\yo 

Similarly, the interference from the simultaneously active D2D-Txs outside the representative 
cluster, xo, at the D2D-Rx of interest can be expressed as: 


'^'I'm\Tx—cluster ^ ^ ^ ^ ll^:' “f ?/|| • (5) 

xS'I’cVaJo 

Denoting the total interference power experienced by the typical device by and recalling 
the serving distance to be r = ||j/o -f xo||, the SIR experienced by the typical device is 


SIR(r) = 


Pdhor- 


Pdhor- 


( 6 ) 


cluster “f -^Tx—cluster 

For notational simplicity, we assume that the system operates in the interference limited regime, 
i.e., the background noise is negligible compared to the interference and is hence ignored. This 
means that the transmit power term cancels in the SIR expression above and can hence be 
ignored, i.e., we can set Pd = 1 without any loss of generality. For a quick reference, the 
notation used in this paper is summarized in Table 


III. Distribution of the Distances 

This is the first main technical section of the paper where we characterize the distributions of 
the distances from the typical device to various intra- and inter-cluster devices. These distance 
distributions will be used in the analysis of coverage probability and ASE in the next section. We 
first focus on the unordered case in which we characterize the distribution of the distances from 
a typical device to a device chosen uniformly at random in a given cluster x G <hc- Using this, 
we characterize the inter-cluster distance distributions for both the content availability cases, and 
serving and intra-cluster distance distributions for the uniform content availability case. We then 
analyze the ordered case in which the distances from a typical device to the devices of a given 
cluster X G <hc are ordered in the increasing order. Using this, we characterize the serving and 
intra-cluster distances to the typical device in the k-closest content availability case. 

Before going into more technical details, we define the functional forms of the probability 
density functions (PDFs) of the Rayleigh and Rician distributed random variables, which will 
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TABLE I: Summary of notation 


Notation 

Description 


Independent PPP modeling the locations of D2D cluster centers; density of $c 


Set of devices inside the cluster centered at a: £ $c 

N = \N^\ 

Total number of devices per cluster (assumed same for each cluster) 

MUM? 

Subsets of denoting the set of possible transmitting and receiving devices 

c Mf, fh 

Set of simultaneously active devices inside the cluster with mean m 


Set of distances from the typical device to the inter-cluster devices 

grro c A/;"'"; m 

Set of simultaneously active devices inside the cluster with mean fh in the cluster Xq £ <I>c 


Set of distances from the typical device to the intra-cluster devices 

'I'm = Elx^(j,^B'^ 

Set of all simultaneously active devices in the D2D cluster network 

II 

Number of maximum possible transmitting devices per cluster, M = N/2 


Scattering variance of the cluster member locations around each cluster center 

Pd 

Transmit power of each device engaged in D2D communications 

a 

Path loss exponent for all wireless links; a > 2 


Channel power gain under Rayleigh fading where hy,^ exp(l) 

Pc;/3;ASE 

Coverage probability; target SIR; area spectral efficiency 


significantly simplify the notation in the rest of this section. 

Definition 1 (Rayleigh distribution). The PDF of the Rayleigh distributed random variable is 
Raypdf(a;a^) =/^(a) = ^exp , a > 0, (7) 

where a is the scale parameter of the distribution. 

Definition 2 (Rician distribution). The PDF of the Rician distributed random variable is 

Ricepdf(a,6;(T^) =/A(a|6) = ^exp (-® ^ 

where /o(.) is the modified Bessel function with order zero and a is the scale parameter. 

A. Distance Distributions in the Unordered Case 

Lets start our discussion with the intra-cluster distances by focusing on the representative 
cluster. Denote by the set {Si}i=i.,M of distances from the typical device to the set of 
possible transmitting devices in the cluster xq E where Si = ||a;o -l-|/|| is the realization 
of Si. The ordering in this case is arbitrary, which means Si E will be interpreted as the 
distance from the typical device to a device chosen uniformly at random from Whenever 
this interpretation is clear, we will drop index i from Si and Si. Characterizing the marginal 
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distribution of S is quite straightforward. Since Xq and y are i.i.d. Gaussian random variables with 
variance a'^, xo + y is also Gaussian with variance 2a‘^. Therefore, S is Rayleigh distributed with 
probability density function (PDF) fsis) = Raypdf (s; 2cr^). However, this does not completely 
characterize because it doesn’t capture the fact that the distances from intra-cluster devices 
to the typical device {||a:o + |/||} are correlated due to the common factor xq. That being said, if 
we condition on the location of the cluster center xq relative to the typical device, the distances in 
the set are i.i.d. since the device locations are i.i.d. around the cluster center by assumption. 
This conditional distribution is characterized in the following Lemma. In the proof, we show that 
instead of conditioning on the location xq, a “weaker” conditioning on the distance uq = ||xo||, 
suffices. Therefore, the statement of the Lemma is presented in terms of vq. 

Lemma 1 (Distribution of i.i.d. sequence Conditioned on the distance uq = ||xo||, the PDF 
of an element S chosen uniformly at random from the i.i.d. sequence S^° is 

fsis\xo) = fsisluo) = Ricepdf(s, i/o;ct). (9) 

Proof: See Appendix ■ 

Remark 1 (Serving and intra-cluster distances for uniform content availability). As discussed 
above, the M elements of are i.i.d. with the distribution characterized by Lemma For 
the uniform content availability case, one of the elements of is chosen uniformly at random 
as the serving distance, and the rest correspond to the distances to the possible intra-cluster 
interfering devices. Since the elements of were “unordered” and the selection of serving 
device was done uniformly at random, all these distances are i.i.d. and follow Rician distribution 
given by Q. The results are stated formally as Corollaries of Lemma below. 

Corollary 1 (Uniform content availability: serving distance). For the uniform content availability 
case, the conditional PDF of the serving distance r = ||xo -f i/o||, conditioned on the distance 
^0 = ||a:o|| between the cluster center and the typical device, is fR{r\i'f) = Ricepdf (r, i^o; cr)- 

Corollary 2 (Uniform content availability: intra-cluster interferer distance). For the uniform 
content availability case, the distances from the intra-cluster interfering devices to the typical 
device {w = ||a:o + y\\,'^y G \ yo} are conditionally i.i.d., conditioned on z/q = ||xo||, with 
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each distance following the PDF given by = Ricepdf (w, oq] a). 

We now look at the distribution of the distances from inter-cluster devices to the typical 
device. Recall that in both the content availability strategies, the inter-cluster interfering devices 
are chosen uniformly at random from the set of transmitting devices fff in each cluster x G <hc- 
Therefore this discussion is applicable to both the strategies. Denoting the distances from inter¬ 
cluster interfering devices of the cluster a: G <hc to the typical device by , it can be shown 
that the elements of are conditionally i.i.d., conditioned on the distance u = ||x|| form the 
typical device to the cluster center x G <hc- h follows on the same lines as Lemma except 
that conditioning here is on i / = ||x|| and not uq = ||xo||. The result is formally stated below. 

Lemma 2 (Inter-cluster interferer distance distribution). Conditioned on the distance v = ||x|| 
between the cluster center x G <hc cind the typical device, the distances from the inter-cluster 
interfering devices to the typical device {u = ||x -1- y\\,'iy G B^} are i.i.d. with each element 
following the PDF given by fu{u\i') = Ricepdf (m, z/; a). 

Proof: Recall that every y G is an independent zero-mean Gaussian random variable. 
Hence conditional on the common distance u = ||x|| between the cluster center x G <hc and the 
typical device, the distances {u = ||x -I- y\\,yy G B^} = are i.i.d. and the conditional PDF 
of each element of can be derived on the same lines as Lemma ■ 

B. Distance Distributions in the Ordered Case 

While the unordered case was sufficient to handle the uniform content availability strategy and 
the inter-cluster interfering distances in both the strategies, we need to consider the “ordered case” 
to handle the serving and intra-cluster interfering distances in the fc-closest content availability 
strategy. For this analysis, consider the set of distances {5'*}*=!:^, denoted by from the 
previous subsection, that denotes the “unordered” distances from the typical device to the set 
of possible transmitting devices in the cluster xq G $c with sampling distribution /s(s|z/o) 
given by Lemma Different from the previous subsection, we order the elements of ° in 
terms of the increasing distance from the typical device and denote them by where 

S{i) < ... < 5'(fc)... < 5 '(m)- In this case the serving distance R corresponds to S(k) whose 
distribution can be derived using order statistics. The result is presented in the following Lemma. 
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Lemma 3 (fc-closest content availability: serving distanee). For the k-closest content availability 
strategy, the conditional distribution of serving distance, conditioned on oq = ||a;o||, is 

fR{rH) = (fc _ “ Fs{r\uo))^-''fs{r\uo), (10) 

with fs{r\u(yj = Ricepdf (r, z/q; O') derived in Lemma^being the conditional PDF of S, and 
Fsi'f'Wo) = 1 — Qi{^, “) being the conditional cumulative distribution function (CDF), where 
Qi(a, 6) is the Marcum Q-function defined as Qfia,b) = te ^/o(at)dt. 


Proof: In Lemma we showed that the unordered sequenee of distanees {Si}i=i-,M is 
eonditionally i.i.d. with eonditional PDF of eaeh element given by /s(r|z/o). The result simply 


follows from the PDF of the order statistie of this sequenee of i.i.d. random variables [32 
eqn. (3)]. ■ 

We next eharaeterize the distanees from the interfering deviees to the typieal deviee in the 
fc-elosest eontent availability strategy. In eomparison to the uniform eontent availability strategy, 
this is more involved sinee the If^ elosest deviee is fixed a priori as the serving deviee and 
henee eannot aet as an interferes To address this issue, we divide the set of simultaneously 
aetive deviees into three subsets, where the serving deviee is loeated at 

a distanee S(k) = ||xo + yo\\ from the typieal deviee, and B(f denote the set of deviees 

that are eloser (farther) to the typieal deviee eompared to the serving deviee. We show that 
eonditional on the distanee uq = ||xo||, the distanees from the typieal deviee to the deviees in 
B(^f are i.i.d. and their distribution is eharaeterized in the Lemma below. The same holds for 
the deviees in This i.i.d. property will play a major role in the exaet analysis of Laplaee 
transform of intra-eluster interferenee in fc-elosest eontent availability ease in the next seetion. 


Lemma 4 (fc-elosest eontent availability: intra-eluster interferer distanee). For the k-closest 
content availability strategy, 

a) the distances from the devices in the set B(f to the typical device, i.e., {tCin = ||a;o + 2/i|}, 
conditionally i.i.d., conditioned on the serving distance r and the distance uq = ||xo|| between 
the cluster center and the typical device, with each distance following the PDF 


fw,S'utin\vo,r) 


/s(«'inko) 

Fs{r\uo) ’ 


0, 


Win < r 
Win > r 


( 11 ) 
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where fsiwinWo) = Ricepdf (wi„, z/q; cr), and Fs{r\uo) = 1 and 

b) the distances from the devices in the set the typical device, i.e., {wout = ll^o + l/i|}> 

are conditionally Ltd., conditioned on the serving distance r and the distance oq = ||xo|| between 
the cluster center and the typical device, with each distance following the PDF 


fw„„,iwont\ao,r) 


/s('»)ont|t^0) 

l-Fs(r|iyo) ’ 

0 


Wont > r 
Wont < r 


( 12 ) 


where fs{wontWo) = Ricepdf (wout, a), and Fs{r\i'Q) = 1 - ^). 


Proof: See Appendix 


IV. Coverage Probability and ase Pereormance 

Using the distanee distributions derived in the previous seetion, we now derive the eoverage 
probability of the typieal deviee and the ASE of the whole network for the two eontent availability 
strategies. We begin our diseussion with the uniform eontent availability strategy. 


A. Uniform Content Availability 

In this subseetion, we foeus on the ease in whieh the eontent of interest for the typieal deviee 
is available at a deviee ehosen uniformly at random in the representative eluster. As evident in 
the sequel, we need the Laplaee transforms of the intra- and inter-eluster interferenee powers as 
the intermediate results for the eoverage and ASE analysis. These are derived next. 

1) Laplace Transform of Interference: We start by deriving exaet expressions and several 
bounds and approximations on the Laplaee transform of intra-eluster interferenee. 


Lemma 5 (Laplaee transform of intra-eluster interferenee). Under uniform content availability, 
the conditional Laplace transform of the intra-cluster interference power given by condi¬ 
tioned on the distance uq = ||xo|| between the cluster center and the typical device is 


M-l 


c 


Itx — c 


(sko) = 


fe =0 


1 -f- sw~°^ 


fw{w\uQ)(lw 




(13) 


with ^ For M ^ m, the above expression reduces to 


-{m — 1) 


sw " 


Iq 1-{- sw 


— fw{w\uo)dw , 


^/Tx-clnster(sko) = CXp 


(14) 
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where fw{'>^Wo) given by Corollary^ 
Proof: See Appendix 


Remark 2. Note that the assumption M 3> m, under which a simplified expression is derived 
in the above Lemma, is applicable when the number of simultaneously active devices per cluster 
is much smaller compared to the cluster size. As discussed in the sequel, this is also the regime 
in which the network performance in terms of ASE will be usually optimized, especially for the 
uniform content availability case. Therefore, the simpler expression will be treated as a proxy 
of the exact expression for the derivation of simpler bounds and approximations. 

While Lemma provides an exaet expression for the Laplaee transform of intra-eluster 
interferenee, it is usually desirable to derive simple but tight approximations and bounds whenever 
possible to draw useful system design insights, whieh we do next. First we derive an approxi¬ 
mation for the Laplaee transform of intra-eluster interferenee under the following assumption. 

Assumption 1 (Uneorrelated intra-eluster distanees). Recall that the distances between intra- 
eluster devices and typical device, denoted by S^°, are correlated due to the common factor 
Xq. However, the coverage analysis can be simplified significantly if this correlation is ig¬ 
nored, which we do as a part of this assumption. More formally, we assume that the serving 
and intra-cluster interferer distances are i.i.d. Rayleigh distributed with marginal distributions 
fniv) = Raypdf (r; 2cr^) and fwi.'vo) = Raypdf (te; 2a‘^), respectively. In other words, we do not 
condition on the distance z/q = H^^oll. which simplifies the coverage analysis by allowing separate 
deconditioning on the serving and intra-cluster interferer distances as discussed in the sequel. 

Under the above assumption, the Laplaee transform of intra-eluster interferenee at the typieal 
deviee is stated as a the following Corollary of Lemma 

Corollary 3 (Approximation). Under Assumption the Laplace transform of intra-cluster 
interference at the typical device is 



( 15 ) 


where fwiw) = Raypdf (te; 2 ct^). 
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By applying Jensen’s inequality on the result of Lemma we also provide a closed form 
lower bound on the Laplace transform of intra-cluster interference in the next Corollary. 


Corollary 4 (Lower bound). The lower bound on Laplace transform of intra-cluster interference 
at the typical device is 


^/Tx-cluste.(s) > exp 


/ (m — 1) 2 27r/a \ 

\ 4cr^ sin(27r/a) / 


( 16 ) 


Proof: See Appendix ■ 

We will use this result along with the lower bound on the Laplace transform of inter-cluster 
interference (will be derived later in this section) to derive a closed form approximation for 
coverage probability and ASE. The results will provide insights into several system design 
guidelines. We now state the exact result for the Laplace transform of inter-cluster interference. 


Lemma 6 (Laplace transform of inter-cluster interference). The Laplace transform of inter¬ 
cluster interference experienced at the typical device, given by is 




^'rn\Tx —cluster 


poo tvt „ 

»=exp(-2xA,y 


1 -I- 


— fu{u\v)du 


kfpfe m 


klrj 


z/dz/j, 
(17) 


with T] = Pqp M ^ m, the above expression reduces to 


^j=0 j\ 


C. 


\Tx —cluster 


» = exp ( - 2xA, p-exp(-m^ ^ ^ 




— fu{u\u)<lv}j^udi?j, (18) 


where fu{u\h>) given by Lemma ^ 
Proof: See Appendix 


Remark 3. Recall that the inter-cluster interfering devices are chosen uniformly at random 
from each cluster in both the content availability strategies, which means that the above result 
is applicable for both the uniform and k-closest content availability strategies. 


Under the assumption M m, a closed form lower bound for the Laplace transform of 
inter-cluster interference can be derived, which is stated next. 

Corollary 5 (Lower bound). Using ( fTS] ) from Lemma the following lower bound on the 
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Laplace transform of inter-cluster interference experienced by a typical device can be derived: 

21 : j a 


C-i (■s) > exp i ^ . 

*”\Tx-ciusterV ) - ^m.{2i: j od)) 


(19) 


sin(27r/a) 

Proof: See Appendix ■ 

2) Coverage Probability: The coverage probability is formally defined as the probability that 
SIR experienced by the typical device exceeds a certain pre-determined threshold /) for successful 
demodulation and decoding at the receiver. It is mathematically expressed as: 


Pc = ER[P{SIR(i?) > (3\R}]. 


( 20 ) 


Using the Laplace transform expressions of intra- and inter-cluster interference powers derived 
so far in this section, an exact expression for Pc is derived in the following Theorem. 


Theorem 1 (Coverage probability). Using Laplace transform of intra- and inter-cluster inter¬ 
ference, given respectively by Lemmas and the coverage probability is 


^oo poo 


Pc = 


JO JO 

where fvoi^o) = Raypdf (z/q; cr^). 


^A,„\Tx-clustec(/^’^“)^/Tx-clu.tec(/^’^“ko)/R(r|z/o)/Vo(^o)drdz/o, (21) 


Proof: From the definition of coverage probability, we have 


Pc = E 


R 


P 


hox^r “ 


>/) 


R 


= Er 


E 


dll'n,\Tx—cluster P ^Tx—cluster 
exp ^ fdv cluster P d'j’x—cluster)) 


R 


where (a) follows from Rayleigh fading assumption. The result now follows from the indepen¬ 
dence of intra-cluster and inter-cluster interference, followed by de-conditioning over R given 
vq using the serving link distance distribution given by Corollary followed by de-conditioning 
over z/q, which is simply a Rayleigh distributed random variable due to the position being sampled 
from a Gaussian distribution in around each cluster center. ■ 

3) Area Spectral Efficiency: The ASE simply denotes the average number of bits transmitted 
per unit time per unit bandwidth per unit area. Assuming that all the D2D-Txs use Gaussian 
codebooks for their transmissions, we can use Shannon’s capacity formula to define ASE = 
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Alog2(l+/9)Pc, where A is the density of the active transmitters and Pc is the coverage probability 
of the typical device. The result specialized to our setup is given in the following theorem. 

Proposition 1. The ASE of the clustered D2D network under uniform content availability is 

ASE = mAclog2(l + /3)Pc, (22) 

where Pc is given by 0 and mXc represents the average density of simultaneously active 
D2D-Txs inside the network. 

Remark 4 (Optimum number of simultaneously active links). Note that there is a clear tradeoff 
between link efficiency and cluster interference. While more active links means potentially higher 
ASE, it also increases interference significantly. ASE can, in principle, be maximized as 

ASE* = max mAc log 2 (l +/3)Pc- (23) 

By solving this ASE optimization problem numerically, we will demonstrate the existence of an 
optimal value of fh that maximizes the ASE in the numerical results section. 

4) Bounds and approximations: After characterizing coverage and ASE exactly, we now focus 
on tight bounds and approximations that will result in easy-to-compute expressions providing 
useful system design guidelines. First, we derive coverage probability under Assumption 

Corollary 6 (Coverage probability approximation). Under Assumption the coverage proba¬ 
bility of the typical device in the uniform content availability case is 

poo ^ 

Pc= / ^/Tx-cluster(/^’^“)>CA„,\Tx-clu.t..(/^’^“)/R(^)dr, (24) 

where fffir) = Raypdf (r; 2ct^), '(^/Tx-ciuster(/^^'^) Laplace transform of intra-cluster inter¬ 
ference under Assumption given by 0. and Laplace transform of 

inter-cluster interference given by lO. 

Proof: The proof follows on the same lines as Theorem except the fact that due to 
Assumption there is no conditioning on z/q (dependence in the intra-cluster distances is 
ignored), which allows to decondition on the serving and intra-cluster interfering distances 
separately. ■ 
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In the numerical result section, we show that this easy-to-compute approximation is also very 
tight. To gain more insights into the behavior of this clustered network, we derive closed form 
approximations for both coverage and ASE in the following corollaries. 

Corollary 7 (Coverage: closed-from approximation). The closed form approximation for the 
coverage probability of a typical device in uniform content availability case is 

Pc —---i— ir~i -■ (25) 

(dvrAccr^m + m - l)/3- + 1 

Proof: The proof follows from the expectation of the lower bounds provided in ( [T^ and ( [T9l ) 
when s = (3r°‘, with respect to the marginal distribution of serving link distance /^(r) = 
Raypdf (r; 2cr^). This is an approximation because of independent de-conditioning over serving 
link and intra-cluster interferer distance distributions. ■ 


Corollary 8 (ASE: closed-form approximation). The ASE in this case can be approximated as 


ASE ~ 


mXr 


(47rAca2m + m - 1)/?: + 1 


(26) 


Proof: The proof follows from the definition of the ASE when Pc is substituted by ( [25] ). ■ 
In the numerical results sections, we show that these closed-form expressions provide sur¬ 
prisingly tight approximations for both coverage and ASE. Based on these two approximations 
(Corollaries [^ and [^, we make the following two observations: (i) ASE and Pc are decreasing 
functions of cr^, which means that the coverage probability and ASE will increase if the devices 
form more dense clusters around cluster centers, as opposed to more spread-out clusters, and 
(ii) increasing density of cluster centers has a conflicting effect on the coverage probability and 
ASE: coverage decreases and ASE increases, which means that the density of cluster centers can 
be increased to enhance ASE as long as the coverage probability remains acceptable. 


B. k-Closest Content Availability 

In this setup, we extend the analysis to the case where content of interest is available at 
closest transmitting device to the typical device. Recall that by tuning the value of k, the content 
can be biased to lie closer (small k) or farther (large k) from the typical device. The best and 
worst case performances can be characterized by setting k = 1 and k = M, respectively. Similar 
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to the previous subseetion, we begin by deriving the Laplaee transform of interferenee. 


1) Laplace transform of interference: Using the distanee distributions derived in Seetion III-B 


the exaet expression for the Laplaee transform of intra-eluster interferenee is derived in the 
Lemma below. As diseussed in the previous seetion, this ease is signifieantly more involved than 
the uniform eontent availability ease due to the way the serving deviee is seleeted. 


Lemma 7 (Laplaee transform of intra-eluster interferenee). Under k-closest content availability, 
the conditional Laplace transform of the intra-cluster interference power given by condi¬ 
tioned on z/Q = llaioll, is = 

M-l Qui 


EE 

n=0 /=0 


, (1 -P) . / , .x 

I I 9mi 9ra “f 1) 


^in(s,r|z/o)Xout(s,r|z/o)’" p 


(27) 


with, 


^in(s,r|z/o) = 


/Cout(s,r|z/o) = 




1 + sw[^^ 
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^ fw^ut (Want I , r) dWout, 

1 + 


where 9 ^ = min(n, k - 1 ), p = pn = 


k-l „ _ , (■ _ y^M-1 


with ^ ^. Q 




n\^ ^i =0 j\ ’ 

is regularized incomplete beta function, fwi^iwifuQ^r) and /wout(^out|z^ou) <^re given by ( |35] ) 
and ( [T^ respectively. Note that here zero to the zero power is defined as one. 

Proof: See Appendix ■ 

Corollary 9 (Best link: 1®* elosest). The Laplace transform of intra-cluster interference for best 
link, i.e., k = 1, under condition offh<^M reduces to 


^iT.-ciustJs,rWo) = exp -(m - 1) 


-00 -a 


Ot \ 

^ , ^^/wout(wzoutkou)dwout , (28) 

1 + SW^nt ) 


where /wout(w^outkou) is given by (fT2l). 


Corollary 10 (Worst link: elosest). The Laplace transform of intra-cluster interference for 

worst link, i.e., k = M, under condition offh^M reduces to 


^/Tx-ciu.te.(suko) = exp -{m - 1) 


f*r —a 

SWin 


1 + SW': 


3^/Win(Win|z^oU)dWi, 


(29) 
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where given by ([^ 


2) Coverage Probability Analysis: Using the serving distanee distribution derived in Lemma[^ 
and the inter- and intra-eluster interferenee Laplace transforms derived respectively by Lemmas 
and 1^ the exact expression for the coverage probability in this case is derived below. 


Theorem 2 (fc-closest content availability: coverage probability). The coverage probability of 
the typical device for the k-closest content availability strategy is 


/ ^/^^\Tx-cluster(/^’^")^/Tx-clu.te.(sU|i^0)/ij(r|l^0)/vb(i^0)drdl/0, (30) 

Jo Jo 

where /R(r|z/o) is the serving distance distribution given by fvoivfj = Raypdf(i/o;a^), and 
the inter- and intra-cluster interference Laplace transforms are given by Lemmas and 


Proof: The proof follows on the same line as Theorem and is hence skipped. ■ 

3} Area Spectral Efficiency: Unlike uniform content availability, where each active receiver 
selects a serving device uniformly at random from its cluster, defining ASE for the fc-closest 
content availability strategy is more involved because different devices may have a different 
value of k for their serving devices depending upon the scheduling strategy. Since we are not 
characterizing scheduling policies explicitly, the information about k for each device is not 
known. Therefore, we derive ASE for a simpler case in which all active receivers connect to the 
k^^ closest device, with k being the same for all the devices. Two special cases of interest are (i) 
k = V. all devices connect to the closest device from the set of possible transmitting devices in 
their clusters, and (ii) k = M: all devices connect to the farthest devices from the set of possible 
transmitting devices in their clusters. The ASE for this setup is given by the following Theorem. 


Proposition 2 (fc-closest content availability: ASE). The ASE for this setup is 


ASE = m\c log2(l + /5)Pc, 


(31) 


where Pc is given by ( [30| ) and mXc represents the average density of simultaneously active 
D2D-Txs. The maximum and minimum ASE is achieved for k = 1 and k = M, respectively. 


Due to the dependence of Pc, given by ( |30l ), on k, ASE is also parameterized by k. As noted 
above, this allows us to characterize the maximum and minimum ASE achievable in a network. 
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4) Bounds and approximations: From Lemma we note that the exact expression of intra¬ 
cluster interference for general k involves two summations, which further complicates the nu¬ 
merical evaluation of the exact coverage probability expression given in Theorem Therefore, 
we propose two simple but tight approximations. In the first approximation, we ignore the effect 
of excluding the closest device (serving) from the field of possible interferers. This allows us 
to assume that the intra-cluster interferers are chosen uniformly at random, resulting in the same 
Laplace transform of intra-cluster interference that we had in the uniform content availability 
case given in Lemma The resulting approximation from coverage probability is given below. 


Corollary 11 (/c-closest content availability: coverage approximation). Using Laplace transform 
of interference in ( [T4| ) and ( [T^ , the coverage probability can be approximated as 


P 

^ c — 


'0 JO 


^Lm\Tx-ciuste.(/^’^“)^/Tx-ciuste.(/^’^“ko)/fl(r|i/o)/yo(^o)drdi/o, 


(32) 


where fR{r\i'fj given by and fvoi^o) = Raypdf (z/q; cr^). 


Proof: The proof follows on the same lines as the proof of Theoremand hence is skipped. 


In addition to the assumption made while deriving Corollary if we also ignore the 
correlation among distances between intra-cluster devices to the typical device (see Assumption 
[^, we can use the simpler result for the Laplace transform of intra-cluster interference given by 
Corollary The resulting coverage probability approximation is stated next. 


Corollary 12 (fc-closest content availability: coverage approximation). Using the results of 
approximation given by and Laplace transform of inter-cluster interference in ( [T^ , the 
coverage probability can be approximated as 


Pr. ^ 


with, fR{r) = 


■^^Tx—cluster Am\Tx—cluster 

Ml ^ 


{(3r°)fR{r)dr 


(33) 


{k-l)\{M-k)\ 
where fsir) = Raypdf (r; 2(T^), and Fs{r) = 1 — exp(—^). 


Fs{rr-\l-Fs{r)r-^fsir), 


In the numerical results section, we will show that both these approximations are remarkably 
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Fig. 2: Coverage probability as a function of number of simultaneously active D2D-Txs in the uniform content 
availability case when ct = 10 and Ac = 150 clusters/km^. 


tight and can in fact be treated as proxies of the exaet result if needed. 

V. Results and Discussion 
A. Uniform Content Availability 

1) Validation of results: We first validate the analytieal results and investigate the tightness 
of various approximations derived for eoverage probability in the uniform eontent availability 
ease. In all simulations, the loeations of devices are drawn from a Poisson cluster process over a 
50 X 50 km^ square region. The eluster eenters are spatially distributed as a PPP, and the deviees 
are normally seattered around them. We set the SIR threshold for the sueeessful demodulation 
and decoding, fi, as 0 dB. Comparing the analytieal and the simulation results in Fig.[^ we note 
that they are a perfeet mateh, thereby validating our analysis. The results also show that the 
approximation for eoverage probability given by Corollary is surprisingly tight. The elosed- 
form approximation given by Corollary also provides a decent lower-bound. 

2) Impact of inter- and intra-cluster interference on the coverage probability: We plot the 
eoverage probability in the presence of (i) only inter-cluster interferenee, (ii) only intra-eluster 
interferenee, and (iii) total interferenee for different seattering varianees in Fig.and for different 
densities of eluster eenters in Fig.|^ It is easy to infer that the eoverage probability of the typieal 
device is strongly dietated by the intra-eluster interferenee. Interestingly, the eoverage probability 
in the absenee of inter-eluster interferenee is independent of the seattering varianee in Fig.[^ This 
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Fig. 3: Coverage probability versus number of simultane¬ 
ously active D2D transmitters for Ac = 150 cluster/km^ 
and different value of a. 



Fig. 4; ASE versus number of simultaneously active D2D 
transmitters for Ac = 150cluster/km^ and different value 
of a. 


is because scattering variance has two counter-balancing effects that cancel each other exactly: 
decreasing scattering variance increases intra-cluster interference by reducing the link distances 
to the interfering devices, and improves the desired link quality by again reducing the distance 
to the serving device. For coverage probability in the absence of the intra-cluster interference, it 
can be seen that the smaller scattering variance provides higher coverage probability, which is 
also true for the coverage probability computed in the presence of both inter- and intra-cluster 
interference. In comparison to the setup of Fig.|^ we decrease the cluster center density and 
increase the scattering variance in Fig.[^ While smaller cluster-center density reduces inter¬ 
cluster interference and hence improves coverage, increasing cluster variance reduces coverage 
(as discussed above in the context of Fig. |^, which has a dominant effect in Fig. 

3) Optimum number of simultaneously active links: In Fig.|^ we present ASE as a function 
of the number of simultaneously active D2D-Txs per cluster. While more simultaneously active 
links may improve ASE, they may also increase interference significantly. Interestingly, it can be 
seen that the optimal value of ASE is the same for a range of scattering variances. This is because 
intra-cluster interference has a more dominant effect on this tradeoff compared to the inter-cluster 
interference. Moreover, as shown in Fig.[^ the coverage probability in the absence of the inter¬ 
cluster interference is independent of the scattering variance which leads to the same optimal 
value for a range of scattering variances. It can also be seen that smaller scattering variances 
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Fig. 5: Coverage probability versus number of active 
D2D transmitters for cr = 50 and different value of Ac- 


Fig. 6; ASE versus number of active D2D transmitters for 
cr = 50 and different value of Ac. 


result in higher ASE, whieh highlights the importanee of short range D2D eommunieation. In 
Fig.[^ we present ASE results for the elustered D2D network with a large seattering varianee for 
the different values of eluster eenter densities. By inereasing the density of the eluster eenters, 
both the density of the simultaneously aetive transmitters inside the whole network, i.e., mAc, and 
the inter-eluster interferenee are inereased. While the former inereases ASE, the latter deereases 
it. However, sinee inter-eluster interferenee doesn’t have a dominant impaet on the eoverage 
probability, and henee on ASE, ASE inereases with the density of eluster eenters. A eommon 
observation in Fig.|^ and Fig.|^ is again the same optimal value for the average number of 
simultaneously aetive D2D transmitters per eluster. As noted above, the intuition behind this 
behavior is that the intra-eluster interferenee, whieh has a signifieant impaet on the tradeoff, is 
independent of both the seattering varianee and the eluster eenter density. 

B. k-closest Content Availability 

1) Best and worst link analyses: After studying the performanee of uniform eontent availabil¬ 
ity strategy in the previous subseetion, we now eompare it with the elosest eontent availability 
strategy, in partieular, its two speeial oases: (i) best link {k = 1) and worst link {k = M) in 
Figs. 0 and We first note that the analytieal and simulation results mateh perfeotly, thereby 
validating the analysis of the fc-olosest eontent availability ease. Next, we notioe that eluster 
















25 



Fig. 7: Coverage probability versus number Of simultane-pig g. simultaneously active D2D 

ously active D2D transmitters for Ac = 150 cluster/km transmitters for Ac = 150 cluster/km^and various values of 
and various values of transmitting cluster size M. transmitting cluster size M. 



- Corollary 11: k=3 

- Corollary 12: k=3 - 

0 simulation: k=3 

- Corollary 11: k=5- 

- Corollary 12: k=5 

^ Simulation: k=5 - 

- Corollary 11: k=7 

- Corollary 12: k=7- 

o Simulation: k=7 


Number of the simultaneously active D2D-Txs 


Fig. 9; Coverage probability versus number of active 
D2D transmitters for n = 10 and M = 40. 



23456789 10 

Number of the simultaneously active D2D-Txs 

Fig. 10; Analytical result for ASE versus number of active 
D2D transmitters for cr = 10 and M = 40. 


size M (number of possible transmitting deviees per cluster) has a conflicting effect on the 
performance of best and worst link cases. This is simply because of “order statistics”: larger 
cluster size M decreases the minimum serving distance (best link) and increases the maximum 
serving distance (worst link). Another interesting observation can be made in Fig.[^ where the 
ASE in the best link case increases (in the range considered) with the number of simultaneously 
active links per cluster, thereby providing “scalability” to the D2D network. 
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2 ) Tightness of the approximations: Recall that the exact expression of intra-cluster interfer¬ 
ence for the fc-closest content availability strategy derived in Lemma [^involves two summations, 
which complicates the numerical evaluation of the exact coverage probability expression of 
Theorem This motivated us to derive two approximations for the coverage probability in 
Corollaries [H] and by approximating the Laplace transform of intra-cluster interference. 
These approximations are plotted along with the exact result in Fig.|^ We observe that both the 
Corollaries provide reasonably tight approximations, which get even better for higher values of 
k. The reason behind this trend is the fact that while deriving these approximations, we ignored 
the effect of excluding the closest device (serving) from the field of possible interferers. The 
effect of this approximation is more prominent at the small values of k. 

3) Optimum number of simultaneously active link: As shown in Figs. and 10, both coverage 
probability and ASE increase significantly when the distance between the typical and serving 
devices is reduced, i.e., the value of k is reduced. More interestingly, we note that the optimum 
number of simultaneously active D2D-Txs that maximizes ASE also increases when the content 
of interest is made available closer to the typical device (smaller value of k). This highlights the 
importance of smarter content placement and scheduling in clustered D2D networks. 


VI. Conclusion 

In this paper, we developed the first comprehensive framework for D2D networks that is 
capable of capturing the fact that devices engaging in D2D may have multiple proximate devices 
any of which can act as a serving device. Modeling the device locations by a Poisson cluster 
process, we considered two content availability strategies that dictate the choice of the serving 
device: (i) uniform content availability: each device chooses its serving device uniformly at 
random from its own cluster, and (ii) k-closest content availability: each device connects to its 
k^^ closest device from its own cluster. Using tools from stochastic geometry, we derived exact 
expressions and several easy-to-use approximations for the coverage probability and area spectral 
efficiency. As a key intermediate step, we characterized the distributions of the distances from 
the typical device to various intra- and inter-cluster devices. To the best of our knowledge, this 
work is also the first to derive these distance distributions for the Thomas cluster process. 

Based on the aforementioned analyses, several interesting system design insights can be drawn. 
The most important amongst them is the existence of the optimal number of D2D transmitters 
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that must be simultaneously activated in each cluster to maximize the area spectral efficiency. 
This can be interpreted as the classical tradeoff between higher interference power and more 
aggressive frequency reuse. Interestingly, we show that this optimal number of D2D transmitters 
is not too sensitive to the system parameters in the uniform content availability case. On the other 
hand, it increases significantly when the content of interest is biased to lie closer to the receivers 
in the fc-closest content availability strategy. This highlights the importance of smart content 
placement and scheduling in these networks. Several other insights regarding the performance 
of this system as a function of cluster center density and cluster variance are also provided. 

This work has several extensions both from the D2D network modeling and cluster point 
processes perspectives. For D2D network modeling, it is of interest to conduct measurement 
campaigns to gather spatial data for various types of user clusters at different locations, such as 
coffee shops, restaurants, airpots, and other hotspot zones, using which we can do data fitting 
to find classes of cluster processes that closely resemble real-world clusters. From the analytical 
perspective, an important extension would be to consider in-band network where devices use 
the same spectrum that is used by the cellular system. This brings forth several interesting 


coexistence issues, especially in the context of cluster point processes p^ . Another interesting 
extension is to incorporate the effects of various content placement and content delivery strategies 
in the proposed model. Incorporating the effect of scheduling in these clustered networks presents 
another avenue of innovation. From the cluster point process perspective, it is important to derive 
the distribution of distances from the typical device to other devices for other point processes of 
interest, such as the Matern cluster process. While the “trends” predicted in this work for D2D 
networks are not likely to change with the choice of the cluster point process, these distance 
distributions will facilitate analysis of more general wireless networks with clustered nodes. 


Appendix 

A. Proof of Lemma 

Denote the location of a device chosen uniformly at random in the representative cluster by 
z = xq + y G and its distance from the typical device by s = ||z|| = ||xo + y|| G M+, 
where xq = {xi,X 2 ) and z = {zi,Z 2 ). Using the fact that y is a zero mean complex Gaussian 
random variable, it is easy to argue that zi = M{xi,a‘^) and 2:2 = A/'(a; 2 ,cr^) conditioned on 



28 


xq = {xi,X 2 ). The conditional joint distribution of Zi and Z 2 can thus be obtained as 

1 / {zi-xiY {Z 2 -X 2 Y' 


/z(2^1,2^2|xo) = 
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Since we are interested in the distribution of S = ||Z||, we define zi = scosO and Z 2 = ssin6', 
where 9 = tan“^ (ff) ■ Using the usual transformation approach, we first evaluate 
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where z/q = y/xl + x^ = ||xo||. Now the conditional marginal distribution of S is fs{s\iyo) = 
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where conditioning on z/q = ||xo|| instead of xq suffices. This completes the derivation of fs{s\iyo). 


B. Proof of Lemma 

The joint distribution of order statistics corresponding to the distances from intra-cluster 
devices which are closer than serving device to the typical i.e., conditioned on S(fc) = r 
and z/Q = ll^oil IS /s(,),...S(,_,)|S(,),.o(s(i),...,S(fe_i)|s(fc),z/o) = , where 


S(i) < S(i+i). Here, the joint distribution is /sn),...,S(fe)|7/o('5(i)7 •••, S(Zc) |i^o) 


fS(fc)(s(fc)|l'o) 


(a) 


M 


M! n /s(s(i)ko)ds(fc+i)...ds(M) 


2 = 1 


k „ . M 

Ml n fs (, S ( i )\ oo ) ■■ Yl /s(^wl^o)ds(fc+i)...ds(M) 

i=l i = k+l 
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k -| 7W n QQ 

= M! J]/g(g(,)|z/o) JJ / /s(s|z/o)ds 

i=l i=fc+l'^^(fe) 

M' ^ 

where (a) follows from the definition of the joint density function of M order statistics [ |34[ 
eqn. (2.10)] with sampling distribution /s(s|z/o), and (b) follows from the symmetry of fsis^i^luo) 

p4| eqn. (2.12)]. Thus, 


(aS)I Ui=ifsis{i)Wo)il - Fs{s(^k)Wo)) 


/5(i),...S(;,_i)|S(fe),i/o(5(l), •••,S(fc_i)|S(fc),Z/o) - 


M-k 


/S(fe)(S(fc)|z2o) 


k-1 


imTT fs{s{i)\i^o) , 


) 'S(i) S 


where (c) follows by using the serving link distribution /s(j^j(s(fc)|z/o) given by ( [TO] ) with r = S(fc). 
Note that there are (fc — 1)! number of possible permutations of the ordered set {s(i)}j=i:fc-i. 
Since the interfering devices are chosen uniformly at random, the permutation term (fc — 1)! does 
not appear in the conditional joint distribution of the unordered set , 


k-1 


/si,...5fc_i|S(fc),!/o(-5l) ■■•5 Sfc-l|S(fc), Z/q) - 


fs(,Si\uo) 

y- Fs{s(k)Wo)' 


^ ^{k) 


(34) 


The product of distributions with the same functional form in the above expression 

implies that the random variables of the unordered set {sj,Vsj < S(fc)} are i.i.d. Therefore, the 
sampling distribution of a set of A; — 1 i.i.d. random variables in the set {wm = ||xo + 2/||} is 


/Wi„(Win|l^0,r) = 


Win <Sik)=r 


0 , 


Win > r 


(35) 


The derivation of fwauti.Wout\k'o,r) follows on the same lines and is hence skipped. 


C. Proof of Lemma [^ 

The Laplace transform of intra-cluster interference at the typical device is 


A^Itx- 


TX —cluster 


[s\oq) =E 


exp(-s hy^Jy + xo\ 

y£B^Q\yo 
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— 


— lUiga^O 


y&B^o\yo 

n 


y&B=^Q\yo 


(b) 


M-1 

E 

k=0 


l + s||?/ + a;o||- 


1 . . ^ (ffi- 

—/Y(l/)d|/ X — 


/r 2 l + s|| 2 / + a:o|| 




p(i<'=fciA'<M-i) 


where (a) follows from the definition of Laplace transform along with the fact that channel gain 
hy^^ ~ exp(l) (Rayleigh fading assumption), (b) follows from the fact that the locations of 
intra-cluster devices, conditioned on xq G are independent and expectation over number of 
interfering devices which are Poisson distributed conditioned on total being less than M — 1, 

where ^ {m-iye ( -( Now under the assumption m -C M, the Laplace transform of 

intra-cluster interference reduces to 


exp ^ - 

(a) 

= exp 


{m — 1) 

— (m — 1) 


l‘ slb + xoll ^ 
R2 1 + sll?/ + a;o||' 


:fY{y)dy 


sw “ 


1 -I- sty““ 


fw(,w\uo)dw), 


(36) 


where (a) follows from the change of variable ||xo -f |/|| w, and converting coordinates from 
Cartesian to polar by using distance distribution given by Corollary 


D. Proof of Corollary ^ 


Using conditional Laplace transform of intra-cluster interference in ( [36| ), the Laplace transform 
of intra-cluster interference can be independently de-conditioned as follows: 

sib + a^oll’ 


"(^Lx-cluster ^ ^ ^ 


(s) = / exp ( — (m — 1) 
/m>2 V 


'R2 1 + 'Sllf/ + ^oll 


—fv(y)dyjfY(xo)dxo 


= I exp ( —(m — 1) 

(b) 

> e 

(c) 


(d) 

= exp 


s z 


1 + s\\z\F^^^^^ ~ ' fYixo)dxo 
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m — 1 2 271 /a 
= exp I —Ti— — —s^ 


sin(27r/a) 

where (a) simply follows from ehange of variable xq + y ^ z, (b) follows from Jensen’s 


inequality, (e) follows from the eonvolution definition, and (d) from Young’s inequality [351. 


E. Proof of Lemma 

The Laplaee transform of the aggregate interferenee from the inter-eluster interferers at the 


typieal deviee, 


^ni\Tx —cluster ' 


(s) = E exp ( - sExe$c\xo KW^ + vW is equal to 


= [ JJ Eb-J JJ E/,^^ [exp(-sVila; + I/ll “)] 

a:e'I>c\a:o 


= E. 


= E. 


[ n Ea-in 


X&^c\xo 

M 

[E( 

k=0 ' 




1 + s||a; + y\\' 


(c) / , 

= exp I — Ac 


(‘^) I o \ 
= exp I — /ttAc 


1 + s\\y + a;| 


M 


F{K = k\K < M) 




A:=0 


1 + s|||/ + X 


v-:fY{y)^y 


klr] 

k ff^k p—fh 


dx 


/ ^ 1 \ k —\ 


k=0 


klr] 


where (a) follows from expeetation over hy^ ~ exp(l), (b) follows from the expeetation over 
number of interfering deviees per eluster, (e) follows from the probability generating funetional 
(PGFL) of PPP [ |3^ , and (d) follows by eonverting from Cartesian to polar eoordinates. Now, 
under the assumption of m M, the Laplaee transform of inter-eluster interferenee is: 




'Pm \Tx—cluster 


(s) = exp ^ — 27rAc / — exp ^ — m 


su " 


'o 1 + su 


—fu{u\iy)du'^^iydi?j. 


F. Proof of Corollary 


The lower bound ean be derived as follows: 




^m\Tx—cluster 


(a) 

[s) > exp ( — 27rA, 


(fe) 


'0 


su “ 


m 


'0 


1 -f SK 


—fu{u\i>)duodo 


su “ 


= exp ( — 27rAc (m / —— —udu) 

Q 1 T su ^ / 


( 38 ) 
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= exp \ 

\ sin( 27 r/a) / 

where (a) follows from the exponential Taylor series expansion and the faet that 1 — exp(—ax) < 
a, a > 0, and (b) follows from the Rieian distribution property that fu{u\u)udu = u, where 
fuiuli') is given by Lemma 


G. Proof of Lemma 

The Laplaee transform of interferenee from disjoint sets of intra-eluster interferers and 

^OUt is >C/Tx-clueter(s>’^l^o) 


= E 




= E 


= E 


exp I —s Vo 11^0+ 2/11 "+ Voll^o + 2/ir 

n exp (-S Vo 11^0 + 2/11““) n exp(-sVol|3^o + 2/ir“) 

1 T-r 1 ^ 


n 


n 


^.0 ^+ ^11^0+ ■+-0 l + s|ko + 2/||- 

yes. “ yeB„“t 


r ffm 

1=0 


n\ ( k — (M — k 


n—l 


lJ\M-lJ \M-lJ Ii_k^{n-gm,gm + l) 


P(i|L<gm) 


/o 1 + sw 


ZI^fWi^{u!inWo,r)dWi, 


1 + SWo “ 


fw.ut{woutWo,r)dwoutj 


n—l-\ 


ICin{s,r\uo) 


ACout/sxIi'o) 


n=0 Z=0 


n!^ 


with g^ = min(n, k — 1), and ^ [m-iye ( -where (a) follows from the faet that 

and B^^^ are two disjoint sets, (b) follows from the definition of Laplaee transform and the 
faet that ~ exp(l), (e) follows from the expeetation over the number of deviees in B^f 
and where the number of deviees in B^f is truneated binomial distribution due to the faet 
that I < min(n, k — 1), along with the faet that distanees from interfering deviees to the typieal 
deviee eonditioned on r and uq in eaeh set are i.i.d., and (d) from the faet that number interfering 
deviees is Poisson distributed with mean rh — 1 eonditioned on the total being less than M — 1. 
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